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Abstract 

We establish various small deviation inequalities for the extremal (soft edge) eigen- 
values in the /5-Hermite and /3-Laguerre ensembles. In both settings, upper bounds on 
the variance of the largest eigenvalue of the anticipated order follow immediately. 

1 Introduction 

In the context of their original discovery, the Tracy- Widom laws describe the fluctuations 
of the limiting largest eigenvalues in the Gaussian Orthogonal, Unitary, and Symplectic 
Ensembles (G{0/U/S}E) [221123] . These are random matrices of real, complex, or quaternion 
Gaussian entries, of mean zero and mean-square one, independent save for the condition that 
the matrix is symmetric (GOE), Hermitian (GUE), or appropriately self-dual (GSE). The 
corresponding Tracy- Widom distribution functions have shape 

i^rvK(i) ~ e^^*' as t -oo, 1 - FtvkI^) ~ e~i^*''" as t ^ cx), (1.1) 

where /? = 1 in the case of GOE, /3 = 2 for GUE, and /5 = 4 for GSE. 

Since that time, it has become understood that the three Tracy- Widom laws arise in a 
wide range of models. First, the assumption of Gaussian entries may be relaxed significantly, 
see [20], [21] for instance. Outside of random matrices, these laws also describe the fluctua- 
tions in the longest increasing subsequence of a random permutation [2], the path weight in 
last passage percolation [11], and the current in simple exclusion [HI [21], among others. 

It is natural to inquire as to the rate of concentration of these various objects about the 
limiting Tracy- Widom laws. Back in the random matrix setting, the limit theorem reads: 
with Amax the largest eigenvalue in the n^n GOE, GUE or GSE, it is the normalized quantity 
^^''^(Amax — '^\fn) which converges to Tracy-Widom. Thus, one would optimally hope for 
estimates of the form: 
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for all n > 1, all e G (0, 1] say, and C a numerical constant. Such are "small deviation" 
inequalities, capturing exactly the finite n scaling and limit distribution shape (compare 
(11. ip ). Taking e beyond 0(1) in the above yields more typical large deviation behavior and 
different (Gaussian) tails (see below). 

As discussed in [TU |T5] , the right-tail inequality for the GUE (as well as for the Laguerre 
Unitary Ensemble, again see below) may be shown to follow from results of Johansson [TT] for 
a more general invariant model related to the geometric distribution that uses large deviation 
asymptotics and sub-additivity arguments. The left-tail inequality for the geometric model 
of Johansson (and thus by some suitable limiting procedure for the GUE and the Laguerre 
Unitary Ensemble) is established in [3] together with convergence of moments using delicate 
Riemann-Hilbert methods. We refer to [T5] for a discussion and the relevant references, 
as well as for similar inequalities in the context of last passage percolation etc. By the 
superposition-decimation procedure of [lO], the GUE bounds apply similarly to the GOE 
(see also [T6]). 

Our purpose here is to present unified proofs of these bounds which apply to all of the 
so-called beta ensembles. These are point-processes on M defined by the n-level joint density: 
for any /? > 0, 

P(Ai, A2, . . . , A„) = ^ n 1^^- - M'e~^'''^ (1.2) 

At /? = 1,2,4 this joint density is shared by the eigenvalues of G{0/U/S}E. Furthermore, 
these three values give rise to exactly solvable models. Specifically, all finite dimensional 
correlation functions may be described explicitly in terms of Hermite polynomials. For this 
reason, the measure (11. 2p has come to be referred to the /9-Hermite ensemble; we will denote 
it by Hp. Importantly, off of /5 = 1, 2, 4, despite considerable efforts (see [9], Chapter 13 for a 
comprehensive review), there appears to be no characterization of the correlation functions 
amenable to asymptotics. Still, Ramirez-Rider- Virag [18] have shown the existence of a 
general (3 Tracy- Widom law, TWp, via the corresponding limit theorem: with self-evident 
notation, 

n^/'(A„.ax(i^/3) - 2v^) TWp. (1.3) 

This result makes essential use of a (tridiagonal) matrix model valid at all beta due to 
Dumitriu-Edelman [5], and proves the conjecture of Edelman-Sutton [6j. As to finite n 
bounds, we have the following. 

Theorem 1. For all n>l, < e <1 and (3 >1: 

^[KrUHp) > 2v^(l + £)) < Ce-^"^'''/^, 

and 
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where C is a numerical constant. 

The restriction to /? > 1 is somewhat artificial. On the other hand, bounds of this type 
cannot remain meaningful all the way down to /5 = 0. Our method in fact applies to all beta 
bounded below, though with the reported C a function of whatever specified minimal beta. 
Keeping (3 > 1 covers the cases of classical interest while allowing for a clearer picture of the 
achieved beta dependence in our estimates (as well as cleaner proofs). 

For completeness we also mention that for e beyond 0(1), the large-deviation right-tail 
inequality takes the form 



For (3 = 1 and 2 this follows from standard net arguments on the corresponding Gaussian 
matrices (see e.g. [E]). For other values of P, crude bounds on the tridiagonal models 
discussed below immediately yield the claim. 

Continuing, those well versed in random matrix theory will know that this style of small 
deviation questions are better motivated in the context of "null" Wishart matrices, given 
their application in multivariate statistics. Also known as the Laguerre Orthogonal or Uni- 
tary Ensembles (L{0/U}E), these are ensembles of type XX* in which X is an n x k matrix 
comprised of i.i.d. real or complex Gaussians. 

By the obvious duality, we may assume here that k > n. When n oo with the n/n 
converging to a finite constant (necessarily larger than one), the appropriately centered and 
scaled largest eigenvalue was shown to converge to the natural Tracy- Widom distribution; 
first by Johansson [H] in the complex {(3 = 2) case, then by Johnstone ^[3] in the real 
{f3 = 1) case. Later, El Karoui [7j proved the same conclusion allowing n/n oo. 

For (3 = 2 and k a fixed multiple of n, a small deviation upper bound at the right-tail (as 
well as the corresponding statement for the minimal eigenvalue in the "soft-edge" scaling) 
was known earlier (see [HIUSI)) extended recently to non-Gaussian matrices in [S]. 

Once again there is a general beta version. Consider a density of the form (11.21) in which 
the Gaussian weight w{X) = e"^'^^/'^ on M is replaced by w{X) = X(/^/'^)('^-"-+'^)+'^ e~^^^^, now 
restricted to M_|_. Here k can be any real number strictly larger than n — 1. It is when k 
is an integer and /? = 1 or 2 that one recovers the eigenvalue law for the real or complex 
Wishart matrices just described. For general k and (3 > the resulting law on positive points 
Ai, . . . , A„ is referred to as the /?-Laguerre ensemble, here for short. 

Using a tridiagonal model for introduced in [5], it is proved in [18]: for k + I > n ^ oo 
with k/u c > 1, 




(1.4) 




(1.5) 
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This covers all previous results for real/complex null Wishart matrices. Comparing (11 .Sp and 
(11. 5p one sees that 0{n^^^e) deviations in the Hermite case should correspond to deviations of 
order (Kny^^{y/K, + y/riy^^e = 0{K^^'^n^^^e) in the Laguerre case. That is, one might expect 
bounds exactly of the form found in Theorem 1 with appearances of n in each exponent 
replaced by n^^^n^^^. What we have is the following. 

Theorem 2. For allK+l>n>l,0<s<l and P > 1: 



p(a_(L,) > (v^+ v^)^(l + .)) < Ce-^^-'''^M^^) 



and 



Again, C is some numerical constant. 

The right-tail inequality is extended to non-Gaussian matrices in j8j . The rather cumber- 
some exponents in Theorem 2 do produce the anticipated decay, though only for e < ^Jnj k. 
For e > ^Jnj K., the right and left-tails become linear and quadratic in e respectively. This 
is to say that the large deviation regime begins at the order 0{^\fnjK) rather than 0(1) as 
in the /3- Hermite case. To understand this, we recall that, normalized by 1/k, the count- 
ing measure of the points is asymptotically supported on the interval with endpoints 
(1 ± \/nj . This statement is precise with convergent ra/zt, and the limiting measure that 
of Mar cenko- Past ur. Either way, \/nj k is identified as the spectral width, in contrast with 
the semi-circle law appearing in the /3-Hermite case which is of width one (after similar nor- 
malization). Of course, in the more usual set-up when c\n < n < (ci > 1 necessarily) 
all this is moot: the exponents above may then be replaced with —f3ne^^'^/C and —(Sn^c^/C 
for e in an 0(1) range with no loss of accuracy. And again, the large deviation tails were 
known in this setting for (3 = 1,2. 

An immediate consequence of the preceding is a finite n (and/or k) bound on the variance 
of Ajnax in line with the known limit theorems. This simple fact had only previously been 
available for GUE and LUE (see the discussion in [15]). 

Corollary 3. Take (3>1. Then, 

Var[A^,.(i/^)] <Cpn-'/\ Var [a^,,(L^)] <CpKn-'/^ (1.6) 

with now constant(s) Cp dependent upon (3. 

The same computation behind Corollary [3] implies that 

limsupn''/^E|Amax(^/3) - 2v^r < oo 
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for any p, and similarly for XmaxiLjs). Hence, we also conclude that all moments of the 
(scaled) maximal Hp and Lp eigenvalues converge to those for the TWp laws (see [3] for 
f3 = 2). 

Finally, there is the matter of whether any of the above upper bounds are tight. We 
answer this in the affirmative in the Hermite setting. 

Theorem 4. There is a numerical constant C so that 

^[XmUHp) > 2v^(l + e)) > c-/'e-^^-''^ 

and 

The first inequality holds for all n > 1,0 < e < 1, and P > I. For the second inequality, the 
range of e must he kept sufficiently small, < e <1/C say. 

Our proof of the right-tail lower bound takes advantage of a certain independence in 
the /3-Hermite tridiagonals not immediately shared by the Laguerre models, but the basic 
strategy also works in the Laguerre case. Contrariwise, our proof of the left-tail lower bound 
uses a fundamentally Gaussian argument that is not available in the Laguerre setting. 

The next section introduces the tridiagonal matrix models and gives an indication of our 
approach. The upper bounds (Theorems 1 and 2, Corollary 3) are proved in Section 3; the 
Hp lower bounds in Section 4. Section 5 considers the analog of the right-tail upper bound for 
the minimal eigenvalue in the /3-Laguerre ensemble, this case holding the potential for some 
novelty granted the existence of a different class of limit theorems (hard edge) depending 
on the limiting ratio n/ n. While our method does produce a bound, the conditions on the 
various parameters are far from optimal. For this reason we relegate the statement, along 
with the proof and further discussion, to a separate section. 



2 Tridiagonals 



The results of [18] identify the general /3 > Tracy- Widom law through a random variational 
principle: 

TWp = snv^-j=j^ f\x)dh{x)-j^ [{nx)f + xf\x)\dxY (2.1) 

in which x ^ b{x) is a standard Brownian motion and L is the space of functions / which 
vanish at the origin and satisfy f'^{x)dx = 1, J^[{f'{x)Y + xf'^{x)]dx < oo. The equality 
here is in law, or you may view (12. ip as the definition of TWp. 
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This variational point of view also guides the proof of the convergence of the centered 
and scaled Amax (of Hp or Lp) to TWp. In particular, given the random tridiagonals which 
we are about to introduce, one always has a characterization of Amax through Raleigh-Ritz. 
In [18], the point is to show this "discrete" variational problem goes over to the continuum 
problem (12.11) in a suitable sense. Furthermore, an analysis of the continuum problem has 
been shown to give sharp estimates on the tails of the (3 Tracy- Widom law (again see |18j). 
Our idea here is therefore retool those arguments for the finite n, or discrete, setting. 

We start with the Hermite case. Let gi,g2, • • •5'n be independent Gaussians with mean 
and variance 2. Let also X2/3, • • • , X(n-i)/3 be independent x random variables of the 
indicated parameter. Then, re-using notation, [5] proves that the n eigenvalues of the random 
tridiagonal matrix 



Hp 



1 

7^ 



gi X{n~l)t3 
X{n-l)f3 g2 X(n-2)/3 



X2/3 gn-1 Xf3 

Xf3 gn 



have joint law (ll.2p j^l Centering appropriately, we define: for t> = (fi, . . . ,Vn) € 



H(v) 



(2.2) 



n-l 



\ 2 

5^ 3kvl + ^ 5^ Xn-kVkVk+l - 2v^ ^ 



k=l 



The problem at hand (Theorem 1) then becomes that of estimating 

Pf sup H{v) > ^/n£] and pf sup H{v) < -y/ne), 



(2.3) 



where we have introduced the usual Euclidean norm ||f||2 = X]fc=i^fc- ^o make the connec- 
tion between Hip) and the continuum form (12.11) even more plain we have the following. 



Lemma 5. For any c > define 



2 

Hc{v) = 3kvl + [Xp{n-k) - E(X/3(n-fc))] V^Vk^ 

V ^ k=l VP fc^i 

n n 

-Cy/nS2ivk+i -Vkf ^^kvl 

A;=0 ^ k=\ 



(2.4) 



^For /3 = 1 or 2 this can be seen by applying Householder transformation to the "full" GOE or GUE 
matrices, and appears to have been used first in a random matrix theory context by Trotter |25j . 
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in which it is understood that Vq = Vn+i = 0. There exist numerical constants, a > b > 0, so 
that 

Haiv) < H{v) < Hb{v) for all v G M", 

granted (3 > I. 

We defer the proof until the end of the section, after a description of the alhed Lp set- 
up. The point of Lemma [5] should be clear. For instance, for an upper bound on the first 
probability in (12. 3p one may replace H by Hf, with any sufficiently small 6 > 0, and so on. 

The model for Ljj is as follows. For k, > n — 1, introduce the random bidiagonal matrix 



B, 



1 



X(3k 

X/3{n-l) X/3(fc-l) 



Xf32 X/3(K-n+2) 

Xl3 X/3(K-n+l) 



with the same definition for the x's and again all variables independent. (The use of x is 
meant to emphasize this independence between the diagonals.) Now [51 shows that it is the 
eigenvalues of = {Bp){BpY' which have the required joint densityo Note that Lp does 
not have independent entries. 
Similar to before, we define 



y/KL{v) 



.T 



-'- \^ 2 2 , -'- \^ ~2 



k=l 



n-fc+l)^fe 



k=2 



+ ^ 5Z Xf3{i^-k+l)Xf3{n-k)VkVk+l - (v^ + 



k=l 



k=l 



The added normalization by ^/K makes for better comparison with the Hermite case. With 
this, and since k > n — 1, to prove Theorem 2 is to establish bounds on the following analogs 
of (E3D: 

P( sup L{v)>y/ne) and P( sup L(v) < —\^e] (2-5) 

^||l,||2 = l ^ ^l|f||2 = l ^ 

Finally, we state the Laguerre version of Lemma [5l (We prove only the latter as they are 
much the same). 



^Once again, a,t f3 — 1,2 this connection had been noted previously (via Householder), see [19] for example. 
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Lemma 6. For c > set 

n-l 



k=l 



k=0 V 



n 



where 



Zk = - - A; + 1)) , Zk = {xl(n^k+i) -P{n-k + 1)) , 

and Yfe = — ^ (x/3(K-fc+i)X/3(n--fc) - E[x/3(K-fc+i)X/3(n-fc)]) ■ (2.6) 
r/ien, /or all [3 > 1 there are constants a > b > so that Laiy) < L{v) < Lb{v) for all 

Proof of Lemma Writing, 

^ n 2 1 

n-l 



n-l 



k=i VP 
shows it is enough to compare, for every w, 

n-l ^ n 



k=l V k=l 



kvl 

k=l ^ k=l 

and 



n-l n-l 

fc+ly 



fc=i vP fc^i VP 

For this, there is the formula 

.,,r(r/2 + l/2) 



r(r/2) 

and, if r > 1, also the bounds 



V - 1 /2 < Exr < Vr. 



(The upper bound is simply Jensen's inequality and holds for all r > 0; the lower bound 
requires more work.) This translates to 

k 1 Ak 

for k < n — 1, j3 > 1. Hence, J{v) < 81{v) for every v. Conversely, if A; < n/2, 
E[X/3(n-fc)/v^ > while ii k > n/2, - l£[Xf3(n-k)/V^ > It follows that 

J{v) > I{v)/16. □ 



3 Upper Bounds 

Theorems [H and [2] are proved, first for the /3-Hermite case with all details present; a sec- 
ond subsection explains the modifications required for the /3-Laguerre case. The proof of 
Corollary [3] appears at the end. 



3.1 Hermite ensembles 

Right-tail. This is the more elaborate of the two. The following is a streamlined version 
of what is needed. 

Proposition 7. Consider the model quadratic form, 

Hb{v,z) = -^^Zkvl-by/^'^ivk+i-VkY -—=^kvl, (3.1) 

k=l k=0 V ^ i^^Q 

for fixed b > and independent mean-zero random variables {zk}k=i,...,n satisfying the uni- 
form tail bound E[e'^^*] < e'^'^^ for all X E ^ and some c > 0. There is a C = C{b,c) so 
that 

Pf sup Hb{v,z) > eV^) < (1 - e-'3/f^)-ie"'5"^'^'/^ 
for all 6 G (0, 1] and n > 1. 

The proof of the above hinges on the following version of integration by parts (as in fact 
does the basic convergence result in [T8]). 

Lemma 8. Let si, S2, ■ ■ ■ , Sk, ■ ■ ■ be real numbers, and set Sk = Yl'e=i ~ 0- further 

ti, . . . ,tn be real numbers, to = ^n+i = 0. Then, for every integer m > 1, 

n 1 " ^ / 1 k+m—l N 

Sktk = — y^^JSk+m-l — Sk-l]tk + ( — [Si — Sk] ] (tfc+l — tfc). 

k=l k=l k=0 ^ e=k ^ 



9 



Proof. For any T^, A; = 0, 1, . . . , n, write 

n n 



fc=l jt=l 



k=0 k=0 
n n 

— S'fc]()f:fc+i — tk) + ^^[Tfc — Tfc_i]tfc. 



fe=0 fc=l 

Conclude by choosing Tk = ^ Y,'e=k''^ Si, k = 0,1, ... ,n. □ 



Proof of Proposition Applying Lemma [8] with = Zk and tfc = f ^ (bearing in mind that 
Vq = Vn+1 = 0, and we are free to set = for A; > n + 1) yields 



n ^ k+m—1 



y^zkvi < —y^\Sk+m-i-Sk-i\vi+y^(— v i^^-^fci) 

k=l k=l k=o e=k 

^ n n 
< —^Am{k - l)vl + ^Am{k)\Vk+l + Vk\\Vk+l - V 



m 

k=l k=0 



where 



Amik) = max \Se — Sk\, for k = 0, . . . ,n. (3.2) 

k+i<e<k+m ^ ' 



Next, by the Cauchy-Schwarz inequality, for every A > 0, 

^ n 1 " " . 1 

^l^^k < -m ^"^(^ ~ + ^ J^K+l - Vkf + Jj-nYl ^UkfiVk+l + VkY. 

Choosing A = h^/n we obtain 

sup Hi,{z, v) < max (^j= A^k - 1) + — ^ [A„(A; - 1)^ + A^k)^] - b^) . (3.3) 

And since whenever (j — l)m + 1 < k < jm, 1 < j < [n/m] + 1, it holds 

A:m{k) V A„(A; - 1) < 2A2„,((j - l)m), 

we may recast (13.31) as in 

sup Hb{z,v) 

\\v\\2=l 

/ 2 , , X 4 , , x2 , (?-l)m+l\ 

< max j= A2m{{j - l)m) + A2„,{{j - l)m) -b^ ' ^ 

i<j<ln/m]+i \my/P ^ ' b^/n(3 ^ ' 
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Continuing requires a tail bound on ls.2m{J) for integer J > 0. By Doob's maximal 
inequality and our assumptions on Zk, for every A > and t > 0, 



Pf max Si>t]< e-^*E [e^^^-l < e-^*+2'=™^' 



Vl<£<2m 

Optimizing in A, and then applying the same reasoning to the sequence —S^ produces 



Pf max \SA >t] < 2e" 

\l<£<2m / 



Hence 



p(A2™(J) >t) <2e-*'/«^"^, (3.4) 

for all integers m > 1 and J > 0, and every t > 0. 
/^From (13 ■4p it follows that 

Pf max A.4(j - l)m) - t^^^' " ^);?+ ^] ) > £v^) (3.5) 

E'"'''"'^^^/ 2 ^ , X 6[(? — l)m + ll eJn 
P ^ A2™( 7 - 1 m) > ^^^^ '— i + ^— 



. '"V^^' ( /3m r 6[(j - l)m + 1] , i2\ 
i=i \ V / 



and similarly 

^ r7=;?^2m J-IM 7ri= )>^^) 3.6 

\l<j<[n/m]+l\b^^np 2^n J 2 J 

< 2 '"y ' cxp i [ ^[^^' ~ + ^] I c 

^— ^ \ ({Arm , /n 



Combined, this reads 



pf sup ifi,(z,w) > (3.7) 

^Ikll2 = l ^ 



/3mne2/l28c _|_ / f \ -l3ben/6Acm 



< \p-prnne / izoc , 

— _ g-/3e6m2/64cy V 1 — 6"'^''^/^"''= / 

which we have recorded in full for later use. In any case, the choice m = [s^^^"^] will now 
produce the claim. □ 

We may now dispense of the proof of Theorem [1] (Right-Tail). Before turning to the 
proof, we remark that if e > 1, one may run through the above argument and simply choose 
m = 1 at the end to produce the classical form of the large deviation inequality (11. 4p known 
previously for /5 = 1,2. 



11 



We turn to the values < £ < 1. The form (12.41) is spht into two pieces, 



Hb{v) = Hb/2{v, g) + Hb/2{v, x), 

Proposition [7] applying to each. 

The first term on the right is precisely of the form (13. ip with each Zk an independent 
mean-zero Gaussian of variance 2, which obviously satisfies the tail assumption with c = 1. 
The second term, Hh/2{v,x)j is a bit different, having noise present through the quantity 
YlkZl{Xp(n~k) - Exi3{n~k))vkVk+i- But Carrying out the integration by parts on tk = VkVk+i 
(and Sk = Xpin-k) ~ ^^/^{n-A:)), wiH producc a bound identical to (13.41) . with an additional 
factor of 2 before each appearance of A2m- Thus, we will be finished granted the following 
bound. 

Lemma 9. For x « X random variable with parameter greater than or equal to one, 

Proof of Lemma This may be viewed as a consequence of the dimension free concentration 
inequalities for norms of Gaussian vectors [13j, but here is an elementary derivation. 

A X of parameter r has density function f{x) = CrX^~^e~^'^^'^ on IR+, requiring us to show 
that 

xPe~^^/'^dx ~ V xPe~^'^/'^dx J 

for any p > 0. The case p = can be done by hand, and we neglect it here. Also, we will 
consider only A > 0, things being quite the same for A < 0. 

Taking logarithms and then differentiating in A, we find that the above inequality (for 
p > 0) is implied by 

/•oo /"OO 

> / / {x-\-y){xyYe-^''-^^''\-y^'''dxdy (3.9) 

/»00 /"OO 

= P / {x^-^y" - xPyP-^)e-'^''-^^''^e-y''^ dxdy. 
Jo Jo 

As for this, let X and Y be positive random variables with density functions c^^A^^'^e"^^"^-*^''^ 
and Cqy'^e~y^^'^ respectively. Now q = p — 1 > —1, and we still have A > 0. It is easy to 
convince oneself that KY < EX, which is exactly the second line of (13. 9p . □ 



Left-Tail. This demonstrates yet another advantage of the variational picture afforded by 
the tridiagonal models. Namely, the bound may be achieved by a suitable choice of test 
vector since 

pf sup Ha{v) < -2C^/^e) < F(Ha{v) < -2C^e \\v\\^ 

\\v\\2 = l ^ ^ ^ 
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for whatever {vk}k=i,...,n on the right hand side. (We have thrown in a constant C for reasons 
that will be clear in a moment.) Simplifying, we write 

^(^Ha{v) < -2Cy/^e\\v\\l^ (3.10) 
< F(^Ha{v,g) < -CV^e\\v\\l^ +P(x(^) < -CV^e\\v\\iy 
where in Ha{v,g) we borrow the notation of Proposition [7| and 

2 

Focus on the first term on the right of fl3.10p . and note that 

F(^Ha{v,g)<-CV^6\\v\\l^ (3.11) 



{g E <) ^ ^v^^ E - E(^'^+i - ^^)' - ^ E 

fc=l fc=l fe=0 ^ *;=! 



with a single standard Gaussian. Our choice of v is motivated as follows. The event in 
question asks for a large eigenvalue (think of y/ne as large for a moment) of an operator 
which mimics negative Laplacian plus potential. The easiest way to accomplish this would 
be for the potential to remain large on a relatively long interval, with a flat eigenvector 
taking advantage. We choose 



k / k 
— A 1-- 

ne V ne 

for which 

n 



k / k \ 

Vk = — All ) ioi k < ne and zero otherwise, (3.12) 

ne V ne/ 



E^^~E^^~^^' '^('"k+i-Vk f —, and ^ A;t;^ ~ n^el (3.13) 



ne 

k=l k=l k=0 k=l 



(Here a ^ h indicates that the ratio a/h is bounded above and below by numerical constants.) 
Substitution into fl3.1ip produces, for choice of C = C{a) large enough inside the probability 
on the left, 

¥[Ha{v,g) < -C^e\\v\\fj < e"^"'"'/^ for ne^^^ > 1. 

The restriction of the range of e stems from the gradient- squared term; it also ensures that 
en > 1 which is required for our test vector to be sensible in the first place. 

Next, as a consequence of Proposition (see fl3.8p ) we have the bound: for c > 0, 

-My.'''] /sE^'^i • (3.14) 

^ k=l ' ' k=l ' 
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With c = C^/ne and v as in (13.121) . this may be further bounded by q-^^'^^^/^ . Here too we 
should assume that rae^/^ > 1. 

Introducing a multiplicative constant of the advertised form extends the above bounds 
to the full range of e in the most obvious way. Replacing e with e/C throughout completes 
the proof. 

3.2 Laguerre ensembles 

Right-Tail. We wish to apply the same ideas from the Hermite case to the Laguerre form 
Lhiy) (for small h). Recall: 

^ n 1 2 

Lbiy) = -=^Zkvl + ^^Zkvl + —=^YkVk+iVk (3.15) 



v^t^ v^fc 

n 1 

-h^/n^^{vk+l -Vkf - h— ^ kvl- 

k=0 ^ k=l 



Here, Zk, Zk and are as defined in (12.61) . and the appropriate versions of the tail conditions 
for these variables (in order to apply Proposition [7]) are contained in the next two lemmas. 

Lemma 10. For x be a x random variable of positive parameter, 

E[e^^'] < e^[>^'l(A+2A2) ^^^^ A < 1/4. 

Proof. With r = E[x'] > and A < i. 



,1 - 2A. 

Now, since x > — | implies log(l + x) > x — x"^, for any A < ;| the right hand side of the 
above is less e''^^^^'^^) as claimed. □ 

Lemma 11. Let x one? x be independent x random variables, each of parameter larger than 
one. Then, for every A G M such that |A| < 1, 

E[e^(xx-E[xx])] < _1_ exp ( 2(1^^2) Mx? + m' + 2AEMEK]]). 

Proof. For |A| < 1, using inequality (13.81) in the x variable, 

/oo 
E[e^[E(x)+.]x]^^(5) 
-00 

where 7 is the standard normal distribution on M. Now, for every s, with (13.81) in the x 
variable, 

]E[gA(E{x)+s)xj < gA(E[x]+s)E[x]+A2[E[x]+s]V2_ 

The result follows by integration over s. □ 



14 



What this means for the present apphcation is that 

E[e^^^], E[e^^^] < e^^' for all real A < (3.16) 

and 

E[e^^'=] < 2e^^^' for all real A with |A| < y/^/2V2. (3.17) 

To proceed, we split Li,{v) into three pieces now, isolating each of the noise components, 
and focus on the bound for sup||j,||2=i Lb/^^v, Z) (the notation indicating (13.151) with only the 
Z noise term present). One must take some care when arriving at the analog of (13. 4p . In 
obtaining an inequality of the form P(Am(J, Z) > t) < Ce^*^/*^ we must be able to apply 
(I3.16P (and (I3.17P when considering the Y noise term) with A = 0{t/m). But, examining 
( 13. 5p and ( 13. 6p shows we only need consider t's of order y/P{y/n + ^jKe)'m. Thus we easily 
get by via 

Pf sup Lfe(t;) > v^e) < Cpf sup Lfe/3(t;, Z) > v^e/c) (3.18) 

\\v\\2 = \ ^ ^l|f||2 = l ' 

< ( ^1 \ -(inme'^IC , ( ^ A -/jy^e/Cm 

with C = C{b), compare (13. 7p . Setting m to be the nearest integer to ^ {^Y^^ puts 
both exponential factors on the same footing, namely on the order of ^ and 

removes all e, k, and n dependence on the first prefactor. Certainly the best decay possible, 
but requires e < ^Jnj k. Otherwise, if e > \Jn/ we simply choose m = 1 in which case the 
second term of (13.180 is the larger and produces decay e~^v^£/c Happily, both estimates 
agree at the common value £ = ^Jnj k. 

Left-Tail. It is enough to produce the bound for P(La(t>, Z) < — Ca/k^I It"] H) for large a, 
given G M" and a C = C(a) as in the Hermite case. Indeed, (13.16p and (13.170 show that 
La(f , Z) and Laiv, y) will follow suit. 
We have the estimate 

v(La{y,Z)<-C^f^e\\v\\^^ (3.19) 

n 

= F(j2i-Zk)vl > ^/p[c^/l^e\\v\\l-a^/^\\Vv\\l-{a/V^)\\Vkv\\l 

k=l 

< exp 

Here we have introduced the shorthand 

n n n 

ll^lll = E^t l|Vt;||HE(^'^+i-^'^)'' llv^^ll2 = E^^^' (3-20) 

k=l k=0 k=l 





\v\ 


\l - a^/n\\'Vv\ 


\l - {a/y/n)\\^/kv\ 


1212 
I2J 


8| 


\v\ 


|4 
U 
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and have also used the fact that fl3.16p apphes just as well to —Z^. In fact, the sign precludes 
any concern over the required choice of A. For the F- noise term, care must be taken on this 
point, but one may check that all is fine given our selection of v below. 

For the small deviation regime, we use a slight modification of the Hermite test vector 
f l3:T2|) . and set 

vt,= ( —k] A fl - —k] with 5 = in/KY''' 
\ne J \ ne J 

for k < ne/5 and ffc = otherwise. This requires e < 5 = [nj kYI"^ in order to be sensible, 
and produces the same appraisals for ||f II2, ||^^||4, ||V||2, and Hv^f II2 as in f l3.13p . with each 
appearance of e replaced by ejb. Substitution into fl3.19p yields 

p(L,(t;, Z) < -Cy^e\\v\ \l) < exp - (/3/8) A^/^^^ _ ^ ^^^3^^)] ■ 



For e > ^JnjK^ notice that the particularly simple choice of a constant v gives 

P(^L<,(1,Z) < < exp(^-(/3/8)Kne2[C- (2a/n) -a]2). 

Combined, these two bounds cover the claimed result, provided that ic'l'^rt^l'^e^ is chosen 
larger than one in the former. Extending this to the full range of e and all remaining 
considerations are the same as in the Hermite setting. 



3.3 Variances 

We provide details for Amax(-f^/3), the Laguerre case is quite the same. (Neither is difficult.) 
Write 

Var [A„,ax(i^,3)] <n \ P (|A^ax(iy/3) - 2v^| > ^e) de\ 
Jo 

and then split the integrand in two according whether Amax < or A^ax > 2y/n. 

First note that our upper bound on the probability that Amax(-f^/3) — < —\pne 

applies to any e = 0(1). Further, from the tridiagonal model we see that Amax stochastically 
dominates {l/^/J3) maxi<k<n dk- Hence, for 5 < we have the cheap estimate P(Amax(-f^/3) < 
—6y/n) < e~^"^^^, and thus 

for all e > 0. This easily produces 
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For the other range, recall that we mentioned at the end of proof for the right-tail upper 
bound that the advertised estimate is easily extended to the large deviation regime (cf. (11.41) ) 
to read 

This results in 

and completes the proof. 

4 (Hermite) Lower Bounds 

Right-Tail. This follows from another appropriate choice of test vector v. To get started, 
write 

PI sup 

V||„||2 = l /V / 

> ^[Ha{v,g)>2^e\\v\\ij^{x{v) < V^e\\v\\iy 

Here, as before, x{v) = (2/\/^) Ylk=iiX/3{n-k) - Hx/3{n-k)])vkVk+i- 

Our choice of v is arrived at by examining the first factor above: with, as in the left-tail 
upper bound, a standard Gaussian q, 



Fi^Ha{v,g) > 2eVE\\v\ 

(2 1^ ri n n s 

{g E S ^ 2v^^ 5^ ^ ay^J2^v,+^ - v.f + ^ ^ /cr;M . 
^ k=l k=l k=0 ^ k=l ^ 

Now the intuition is that the eigenvalue (of a discretized —d?'/dx'^+ potential) is being forced 
large positive, so the potential should localize with the eigenvector following suit. 
Let then 

Vk = V^k A (l — \^k) for k < e^^^"^ and otherwise 0, 
where we will assume that n > e^^/^ > e^^^"^. With these choices we have 



\v 



|2~||w||4 ^, \\Vv\\l ^/e, ||v^f||2~-, 

^y e e 



(recall the notation from (13.201) ) and thus the existence of a constant C = C{a) so that 
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Similarly, returning to the second factor on the right hand side of (14. ip and invoking the 
estimate (13.141) we also have 



2^ < -/3n.3/Vc 



for the same choice of v. And granted ne^/^ > 1, it follows that P(x(f) < v^^l 1"^! I2) — 
1 — e^^/*" throughout this regime. That is, 

Pf sup H{v) > ^/^e) > whenever ne^^^ > 1. 

When n5^/^ < 1, write 

P( sup H{v) >V^6)>f( sup H{v) > V^eo) > ^ > 77^6-^^"^''', 
V||t,||2=i / ^||t'||2=i ^ Oe'^'^ oe'^^ 

where £0 = n~'^^'^ < 1 to produce the advertised form of the bound for all n and e. 

Left-Tail. This relies heavily on the right-tail upper bound. The first step is to reduce to 
a Gaussian setting via independence: for whatever 6 > 0, 

pf sup H2b{v) < -^/^e) > pf sup Hb{v,g) < -2y/^e) pf sup Hb{v,x) < Vne). 

\\v\\2 = l ^ \\v\\2 = l ^ \\v\\2 = l ^ 

Here we also use the notation of the proof of Theorem [T] (right-tail), from which we know 
that 

pf sup Hb{v,x) > Vne) < Ce"""'^'/^. 

\\v\\2 = l ' 

(As /3 > 1 we are simply dropping it from the exponent on the right at this stage.) Hence, 
if as we regularly have start with an assumption like ne^l"^ > C*^ > 1, it follows that 



P 



f sup E2h{v) < -y/^e) > (1 - e-^)pf sup Hbiv,g) < -2y/^e). 

\\v\\2 = l ' \\v\\2 = \ ' 



Turning to Hb{v,g) we make yet another decomposition of the noise term. Let L be an 

L 



integer (1 < L < ra) to be specified. Set Sl = J2k=i 9k, and 



Vk = gk- j^Sl for 1 < k < L, r]k = Qk ioi L < k < n. 

Note that the family {rik}k=i,...,n is independent of 5"^. If the procedure of Proposition [7] could 
be applied to Hbiy, rj), we would have an event of probability larger than 1— Ce""^^''^/*" (again 
we simply drop the beta dependence at this intermediate stage) on which 

^ n n ^ n n 

k=l k=l k=l k=l 



Since we are still working under the condition ne^/^ > C^, this is to say that there is an 
event of probability a least 1 — 1/e, depending only of the rj^s, and on which 

^ L h ^ " 

k=i k=i k=i 

for every G M". If we now choose L + 1 > Qne/b, we have further 

^ L L n 

k=l k=l k=L+l 

on that same event. Note this choice requires e < 6/6; it is here that the range of valid 
epsilon gets cut down in our final statement. In any case, putting the last remarks together 
we have proved that 

Pf sup Hbiv,g) < -2^/^6) > (1 - e"^)P (Sl < -3^/nPe) 

\\v\\2 = l / \ J 

and so also 



pf sup H2b{v) < -V^e) > ^ 

\\v\\2=l ^ ^ 



again under the constrains ne^/^ > and e < b/Q. The last inequality follows as Sl is a 
mean-zero Gaussian with variance of order {ne)~^. 
The range ne^/^ < is handled as before, 



P(^ sup^iJ26(^) < -V^e) > P( snp^H2biv) < -v^£o) > > ^ 



where Eq = (C^/n)^/^. As Eq must lie under 6/6, this last selection requires n 
but smaller values of n can now be covered by adjusting the constant. 

It remains to go back and verify that P(sup||^||2=i Hf,{v, rj) > ^Jne) < Ce""^^''^/*" . The only 
reason that Proposition [7| cannot be followed verbatim is that the r/^'s are not independent, 
the first L of them being tied together through Sl. We need the appropriate Gaussian tail 
inequality for the variables 

I 

j=k 

and, comparing with (13. 4p . shows that an estimate of type P(Am(A;,r7) > t) < Ce~*^/*""^ 
suffices. But 

^'no = ^9o^{ir\L-kN L)Sl, 

j=k j=k 

and so 

¥{A^{k,r]) >t) < ¥{A^{k,g) > t/2) +¥{mSL > t/2). 



Arri{k,r]) = max 

k<e<k+m 
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The first term we have aheady seen to be of the required order, and the second is less than 
g-Lt^/sm^^ Since we only apply this bound in the present setting when L = Cne > Ce~^l'^ 
and m = (the choice made in Proposition [7]), we have that ^{raSi^ > t/2) < e"*^/*^™-, 

and the proof is complete. 

5 Minimal Laguerre Eigenvalue 

While not detailed there, the results of [18] will imply that 




whenever K,n — >■ oo,K/n — * c > 1. This appraisal was long understood for the minimal 
eigenvalue of L{0/U}E, and has recently been extended to non-Gaussian versions of those 
ensembles in [8]. The condition n/n —>■ c > 1 keeps the limiting spectral density supported 
away from the origin, resulting in the same soft-edge behavior that one has for Amax- If 
instead n — n remains fixed in the limit, one has a different scaling and different limit law(s) 
for Amin, the so-called hard-edge distributions. Granted the existence of the "hard-to-soft 
transition" for all /5 > (see [1] and [17]) it is believed that (15.11) holds as long as k — n ^ oo, 
but (to the best of our knowledge) this has not been explicitly worked out in any setting. 

We only consider the analogue of the right-tail upper bound for Amin and have the fol- 
lowing. 

Theorem 12. Let (3 > 1 and n > n + 1. Then, 

p(A™„(L,) <{V^- V^)\1 - e)) < Ce~^i^n)^"^V-^-^)e^'^/c^ (5_2) 

for a numerical constant C and allO < e < v^(tt^'^ A a^n"^/^) in which a = 1 — y^n/n. 

According to (15.11) . the deviations are of the order of {^/Kny^^(^/K — y/nY^^e, which ex- 
plains the exponent in (15.21) . Our condition on e is certainly not very satisfactory, although 
still sensible to the fluctuations in (15. ip . One would hope for the range of e to be understand- 
able in terms of the soft/hard edge picture — what we have here arises from technicalities. 
On the other hand, if we place an additional, "soft-edge" type, restriction on k, and n, we 
obtain a more natural looking estimate. 

Corollary 13. Again take (3 > 1, but now assume that k, > en for c > 1. The right hand 
side of Ii5.2\) may then be replaced by Ce~^^'^^''^/^ for a C = C{c), with the resulting bound 
valid for all < e < 1. 
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The last statement should be compared with Corollary V.2.1(b) of [8], which applies to 
classes of non-Gaussian matrices. 

As to the proof, we proceed in a by now familiar way. We first set 



^/kL{v) = v'^ ({-/k - \fnf - Lf^v. 
Then, after a rescaling of e, we will prove the equivalent 

Pf sup L{v) > a''/^^/^e) < Ce-'^"^'^'/^ for e < mm{a^^/'',a^/''n-^/'). 



Similar to the strategy employed above, a series of algebraic manipulations shows that we 
can work instead with the simplified quadratic form 

^ n 1 " 2 



n— 1 2 



k=l 



(The condition k > n + 1 in Theorem [12] is used in passing from L to L' .) 

We remark that under the added condition n > cn for c > 1, a is bounded uni- 
formly from below and 'E\xi3{K-k+i)Xi3(n~k)\ is bounded below by a constant multiple 
of \/n — k. Hence, the deterministic part of L' is bounded above by a small negative mul- 
tiple of X]fc=i ("^fe+i + + X]fc=i ^"^fc- The proof of Corollary [T3] is then identical to 
that of the right-tail upper bound for Amax(-^/3)- 

Back to Theorem [12] and a's unbounded from below, we begin by rewriting the noise 
term in L' as times 

n n n—1 

Y,i-Zk)vi + Y.{-Zk)vi + 2Y,i-yk)vkVk+i 

k=l k=2 k=l 

n n n—1 



in which 



k=l k=2 k=l 



(X/3(K--fe+l) - X/3(n-/fc))^ - - X/3(n-fc))^] 



(with the convention that Xo = 0)- The idea is the following. For moderate k, Var(f/fc) = 
0(q;^), and thus it is this contribution to the noise which balances the drift term X]fc=i ^^l- 
Also, one may check that in the continuum limit the optimal v is such that |ffc-|-ffc+i| = o(l), 
and so the Z and Y terms should "wash out" . 
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We complete the argument in two steps. In step one, we simply drop the Z and Y terms 
and apply the method in Proposition 7 to the further simplified form 



Tl TT\ TT \ 2 \^ ^[Xl3{K,-k+l)Xp{n-k)\ , , x2 a ^ 2 

^(-' f^) = 7^ - E (-'^^i 7^ E 

— 1 k — 1 /l — 1 



(5.4) 



Even here we loose a fair bit in our estimates (resulting in non-optimal on e) due to the 
variable coefficient in the energy term. Step two shows that, under yet additional restrictions 
on e, the Z and Y noise terms may be absorbed into L{v, U). 

Step 1. We wish to prove P(sup||^||^;^ L(f , ?7) > a'^/^^/ne) < (7e-/3"£^''Vc some range 
of e > 0. (The optimal range being < £ < [K/ny/'^a^^^ .) A first ingredient is a tail bound 
on the Uk variables, for which we first bring in the following. 

Lemma 14. (Aida, Masuda, Shigekawa fl^) Given a measure t] on the line which satisfies 
a logarithmic Sobolev inequality with constant C > 0, there is the estimate 

r A(F^-E[F^])^^ < 2e«^^'^[^]' whenever \\\ < 
J 16C 

for any 1-Lipschitz function F . 

As a consequence, we have that: 

Corollary 15. Let x o-nd x be independent x random variables (each of parameter larger 
than one) and set U = {x~ xY ^'^'^ ^ — "^[x " x\- There exists a numerical constant C > 
such that 

for all real A G (-1/C, 1/C). 

Indeed, by the general theory (see Thm. 5.2 of [13] for example) the distribution of the 
pair (x, x) oii ^ satisfies a logarithmic Sobolev inequality. The lemma then applies 
with -F(x, y) = X — y. In our setting, we record this bound as 



E[e^^^] < Ce^"^^' for |A| < y^/C and al = E[U, 



Picking up the thread of Proposition 7, the variable coefficient in the energy term of 
L(f , U) is dealt with by applying the Cauchy-Schwarz argument with A = defined by 

A, = ^!^^^^^^±#^, A: = l,...,n-1, (5.5) 



compare (3.3). Schematically, we are left to bound 

In, I m.l 



V P( 7=A„(jm, U) V — — A„(jm, Uf > ^jm + ea'/^V^) (5.6) 
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for our choice of integer m. The Am(-, U) notation stands in analogy to that used in Section 
3. Note we have taken the hberty to drop various constants and shifts of indices in the above 
display (which are irrelevant to the upshot). 

Here the dependence of and on the relationship between n and k comes into 
play. While at the top of the form everything works as anticipated, these quantities behave 
unfavorably for k near n. For this reason we deal with the sum fl5.6p by dividing the range 
into j < n/2m and j > n/2m with the help of the appraisals: 

l<k< n/2, ^ ^ f v^/C, 1 < A; < n/2, 
n/2<k<n. ~ [ gV^^T^, n/2 <k <n. 

Restricted to j < n/2m (and hence substituting crj^ = Ca^, \jm = Vn/C), the sum 
(15. 6p can be bounded by Ce~^'^^'^'^/^ upon choosing m = [e^^/'^a^'^/^]. This holds for all 
values of e so long as the choice of m is sensible, requiring that e > a~^/^n~^. But this is 
ensured if k > n + 1 and e'^/'^n > 1 (the former having been built into the hypotheses and 
the latter we may always assume). 

On the range j > n/2m the e term on the right hand side within the probabilities is of 




no help, and we use, along with cxj^ < C and Xjm < a^n — jm/C, the crude estimates 

n/2m<j<n/m j>n/2m 



and 



J2 P(A^(jm,f/)2>/?«2^A,„^) < C 2^ e 



n/2m<j<n/m ^ l<j<n/2m 



< C(l + (m/na6))g-/3a3(„/^)i/2/c_ 

The choice of m = being fixed, we can bound each of the above by the desired 

Q^-Pne^^^/c Qj^iy ]jy restricting e to be sufficiently small. The first estimate requires e < a'^^^^, 
the second requires in addition that e < a^^^n"'^^^ (and again uses ne^^'^ ^ !)• 
In summary 



P 



( sup L{v, U) > «^/=^v^£ ) < Ce-^"^'''/^ if < £ < min(a20/3, a^/s^-VS), (5,^ 



It is perhaps worth mentioning here that the bounds on and for the range k > n/2 
introduced in (15.71) may be improved slightly, though not apparently with great effect on the 
final result. 

Step 2. To absorb the Z,Y noise terms, we show that L'{v) < L{v,U) + S{Z,Y,v) 
with a new form L[v, U) comparable to L[v, U), and an "error" term S for which we have 
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F{S > a^/^y^e) < Ce-^"^'^'/^, at least for some range of e > 0. What follows could almost 
certainly be improved upon. 
Define, for k = 1, . . . ,n — 1: 

ttk = —Xk for k < a^n, = — for k > a^n. 
4 16 y/n 

(Recall the definition of from (15.51) .) Then, an application of the Cauchy-Schwarz in- 
equality yields: for all v of length one, 

n-l on 7^2 



^ ^ k=2 k=l k=l 



13 ak 



A similar estimate applies to Yl^lJ^^YkVkivk+i —v^). Accordingly, 

L'(v) <L(v,U)+ max + max — ^ 

l<k<n-l jjak l<fc<n-l jjak 

with 

n n— 1 2 

A:=2 k=l ^ k=l 

Obviously, the arguments of step 1 apply to L{y,U). 

Finally, with either Z^^i or Y^, Lemmas 10 and 11 imply that 



max ^>ea4/3^^ <c<y 

\l<k<n-l jSttk J ^ 



n 



provided say e < 1. Since it may be assumed that a < 1/2 (otherwise we are in the easy 
regime covered by Corollary 13), we have the bound = < \fnjC for k < a^n < n/2 
and so also 

k=l 

by considering the sums over k < a^n and k > a^n separately. If now e < a^'^^^ (still 
keeping in mind that e^^'^n > 1), the right hand side is less than Ce~^"''^'^^^^'-^ . Adding this 
new constraint on e to those stated in (15. 8p completes the proof. 
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